Persistence in a Simple Model for the Earth's 
Atmosphere Temperature Fluctuations 

Arturo Berrones 

Posgrado en Ingenieria de Sistemas 
Universidad Autonoma de Nuevo Leon 
AP 126 - F, Cd. Universitaria 
San Nicolas de los Garza, NL 66450, Mexico. 

February 2, 2008 

Keywords: climate models; stochastic processes; long memory effects; long- 
range correlations; complex systems. 

Abstract 

The effect caused by the presence of a number of distinct time 
scales in a simple stochastic model for the Earth's atmosphere temper- 
ature fluctuations is studied. The model is described by a dissipative 
dynamics consisting of a set of coupled stochastic evolution equations. 
The system shows features that resemble recent observations. In con- 
trast to other approaches, like autoregressive models, the fluctuations 
of the atmosphere's temperature depend on parameters with clear 
physical meaning. A reduced version of the model is constructed and 
its temporal autocovariance function is explicitly written. 

1 Introduction 

Interaction among processes with several length and time scales is common 
to a variety of complex systems. For instance, the long-range temporal cor- 
relations found in signals from a variety of fields can be associated with an 
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interplay of a number of time scales [H |2l El IH E] • In particular, it is an ex- 
tended belief that the persistence observed in the temperature fluctuations of 
the Earth's atmosphere is a consequence of its feedback with slower dynam- 
ical components in the climate system, like the oceans and Earth's surface 
Persistence at short time scales is related to the everyday life observation 
where similar weather conditions are likely to be experienced over a given 
region in a time scale of a few days. The existence of these short-term corre- 
lations make weather forecasting possible. The climate's persistence is also 
found for larger time scales, however its characterization is a more difficult 
task [B]. In some recent experiments, the temperature records from different 
places around the globe have been analyzed. The observations indicate the 
existence of universal power-laws C{t) ~ describing the correlations of 
the temperature fluctuations around its mean seasonal value. Although there 
is some disagreement about the value of the exponent r, it has been firmly 
established that the persistence in temperature fiuctuations can indeed be 
characterized by power-law autocorrelation functions [0]. It has been re- 
ported by some authors that for time scales that range from ~ 1 to ~ 25 
years, correlations measured on data from meteorologic stations placed on 
small islands decay with an exponent r ~ 0.4, while for continental stations 
data are closer to r ~ 0.7 [7j. According to other authors :8], the persistence 
is even more pronounced for the oceanic regions (being roughly characterized 
by a 1/f noise), while in the inner continents r ~ 1 (in terms of the power 
spectrum, a white noise at low frequencies), with a transition region in the 
coastal zones in which r ~ 0.7. 

The emergence of long-range temporal correlations is a non-trivial feature 
that can be used to test models of the Earth's climate [7j. The understanding 
of the long-range temporal correlations is fundamental because they charac- 
terize the interaction among the different climate components [S]. There is 
some controversy with respect to the description made by large scale models 
of the atmospheric temperature variability |Hl IHl UHl HH E] • The purpose 
of the present Letter is to introduce a conceptual stochastic model for the 
fluctuations of the Earth's atmospheric temperature, or more precisely, its 
radiated energy (for a survey on conceptual climate models see, for instance, 
Imkeller and Monahan [121 )• ^^^^ shown that the model displays qual- 
itative features that closely resemble observations. In contrast with other 
simple stochastic models of atmospheric temperature fluctuations, like au- 
toregressive models [1311111, all the parameters of the presented model have 
direct physical interpretation. Additionally, as it will be discussed below, the 
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model introduced here also reveals statistical features that are reminiscent of 
recent observations on the spatial structure of the climate network. There- 
fore, this work is intended to be a contribution towards the construction of 
realistic and unexpensive algorithms for Earth's climate simulation. 

The model is based on energy balance On Earth, as in other plan- 

ets with a solid crust, the influx of solar radiation is balanced by the outflow 
from the surface and the atmosphere. In the simplest description of this 
process, the Earth is treated as a single point. Let us denote by y and x 
the global averages of the radiation emitted by the atmosphere and by the 
surface (oceans and land), respectively. A fraction of the Sun's total radi- 
ation is immediately reflected back into space and another is absorbed by 
the atmosphere. The remainder of the flux is transmitted through the atmo- 
sphere and reaches the surface, which in turn absorbs some of the radiation 
and reflects the rest. The radiation absorbed by the surface is then radiated 
back as heat. The surface radiates the absorbed energy in the infrared (IR) 
region of the spectra. It turns out that the atmosphere is not transparent 
to IR radiation, essentialy due to the presence of the so-called greenhouse 
gases. Let us denote by a the fraction of IR radiation absorbed by the at- 
mosphere. A fraction h of the total radiation absorbed by the atmosphere is 
directed towards the surface and the rest is finally lost into space. All these 
considerations are put together in the well known zero dimensional energy 
balance model: 

y = A + ax, (1) 

x = B + by, 

where constants A and B are the net contribution made to y and x by 
the solar radiation flux, taking into account that some heat is removed by 
water evaporation from the surface. The constant a is called the infrared 
(IR) absorption coefficient. All the constants in Eqs. (0) are calculated 
by averaging over a year and over the entire Earth's surface. Assuming a 
blackbody process, the average atmosphere's temperature is given hj y = 
uT"^, where u is the Stefan-Boltzmann constant. In spite of its simplicity, 
the zero dimensional energy balance model is capable of predicting with very 
good accuracy the mean Earth's surface temperature. Another interesting 
prediction of the zero dimensional energy balance model is the increment of 
the mean temperature as the coupling parameters a and b grow. 

In this Letter a spatially extended and time-dependent generalization of 
model (HI) is introduced. The Letter is organized as follows: in Sec. |2] the 
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model is introduced and the temporal and spatial correlations are discussed 
numerically. A reduced version of the model is constructed and formally 
solved in the framework of the Langevin approach. A discussion of the sta- 
tistical properties of the solution is given. Some conclusions and future di- 
rections are discussed in Sec. El 

2 The Model 

Model is generalized by the assumption that energy balance is satisfied 
locally and a transient time is necessary in order to achieve a stationary 
state. A set of cells is considered. In each cell, atmosphere interacts 
with the surface through the local atmosphere's IR absorption coefficient 
a„ and the local fraction of heat that the atmosphere returns towards the 
surface, Each component, atmosphere and surface, has its own intrinsic 
local response time. The cells of each component interchange radiation via 
a diffusive process. The model is written as 

ijn = diAyn - XnlVn " (^n + a^Xn)] + £„(t), (2) 
Xn = d2AXn - Jn[Xn " (5„ + bnVn)]- 

In this equation (A„)^^ and (7n)^^ are the local response times of the atmo- 
sphere and the surface, respectively. The symbol A is the discrete Laplacian 
and di, d2 represent the diffusion coefficients of each component. The term 
Snit) is a Gaussian white noise, without correlations between different cells. 
The meaning of the rest of the terms follow from the zero dimensional en- 
ergy balance model (P). In particular, i/n gives the radiation emitted by the 
atmosphere in the site n at a given time. The atmosphere is expected to 
have shorter intrinsic response times than that of the surface. The radiation 
emitted by the surface at time t in cell n is represented by x„. The noise 
reflects the more rapid variations, or weather. Periodic boundary conditions 
are taken. The constants A„, 7„, An, Bn, an and bn are assumed to be inde- 
pendent variables, such that averages over index (in the limit A^ oo) give 
the corresponding values of the parameters of the zero dimensional energy 
balance model. A, B, a, b; and the effective inverse response times for the 
atmosphere and the surface. The flrst important thing to notice with model 
Q is that it recovers the zero dimensional energy balance model. This can be 
seen by averaging Eq. Q over the cells in the limit N ^ oo. The following 
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reduced version of Eq. Q is obtained: 

y = -X[y-{A + ax)]+eit), (3) 
X = — 7[x — {B + by)]. 

From the fact that Eq. Q represents an overdamped dynamics in a parabohc 
potential, the system converges to a stationary state with mean value given 
by the solution of the zero dimensional energy balance model. This result 
shows that model Q is capable of representing with good approximation the 
mean behavior of energy in the coupled atmosphere - surface system. Now 
it will be argued that Eq. (0) can give realistic descriptions of temperature 
fluctuations as well. A large variability over the intrinsic characteristic times 
of the surface around the globe is expected to exist, as a consequence of the 
different response times present in the geosphere, hydrosphere, cryosphere 
and biosphere. At first instance, this situation is modeled by treating the 
7„'s like independent random variables taken from a uniform probability 
distribution. On the other hand, under the basis of the relative homogeneity 
of the atmosphere composition, it will be assumed by now that the response 
time to perturbations of the atmosphere is the same in all cells. In what 
follows, the value A„ = A = 1 is used. Under these assumptions, the time 
unit is defined as the mean atmospheric response time. The model represents 
the interaction between a hierarchy of time scales, ranging from minutes to 
days to geological times. From this point of view, it would be therefore 
reasonable to assume that 1/ A lies in the intermediate scales, roughly in the 
range from weeks to months. This and other important aspects about the 
definition of the parameters in the model are intended to be refined by the 
author in the near future by close cooperation with climate experts. At this 
point the main goal is limited to explore the capabilities of the model 0, in 
order to give qualitatively realistic statistical descriptions of the temperature 
fluctuations present in the Earth's atmosphere. 

Figure ^a) is a log-log plot of the power spectrum of the time series 
of the atmosphere's temperature in a particular cell. This time series is 
obtained from the corresponding time series for under the assumption 
of a blackbody process. The parameters a„, bn, An and are left constant 
and set to their experimental global averages: a„ = a = 0.96, bn = b = 0.61, 
An = A = 179.36 W/m2 and Bn = B = 47.82 W/m^ [15]. The values 7„ are 
drawn from a uniform distribution on the range (0,0.2). The noise values 
en{t) are uniformly distributed over the interval (—1, l)W/'m'^. The diffusion 
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coeffiecients are taken as di = d2 = 1- System size is = 50. The power 
spectrum is consistent with a power law at low frequencies, with exponent 
~ —0.25. This exponent of the power spectrum implies a power-law decay 
of correlations at large times, C{t) oc t~^, with r ~ 0.75. An alternative 
approach is given in Fig. ^b). In order to characterize the correlations, the 
sum of the temperature values on time at a given cell is studied, resulting in: 



where T„(r) stands for the atmosphere's temperature at cell n and time r. 
The signal Yn{t) is then compared to a random walk. In Fig. ^b) the 
standard deviation of Yn{t) as a function of time is plotted. For large times 
(roughly greater than 50 time units), a ~ t^-^^^, which implies a power-law 
decay of the autocorrelation function at large times, with exponent ~ 0.7. 
This result is consistent with Fig. ^^a). 

An analysis of the radiation spectrum has been carried out for the same 
experimental setup as above. Results indicate that radiation and temper- 
ature spectra are basically equivalent, differing at most by a normalization 
factor. This numerical finding is consistent with recent studies that suggest 
that a strongly correlated signal preserves its correlation properties after an 
even polynomial transformation jl7j . 

The model is also capable of showing scaling behavior for surface tempera- 
ture. As an example, the power spectrum of surface temperature fluctuations 
of an induvidual cell is shown in Fig. |21 The system has the same parameter 
values as before, but d2 = 0. Notice that with this choice of parameters, the 
coupling with the atmosphere is essential for the emergence of scaling in the 
surface temperature. 

The model displays an interesting spatial structure. In Fig. |2] a case is 
considered in which the system size is = 200 and the other parameters 
are the same as in the case presented in Fig. ^ The spatial autocorrelation 
function C (n) is inferred from the power spectrum of the temperature values 
vector at a fixed time. A least-square fit of the power spectrum to a power- 
law function indicates that S{k) ~ k^^'^"^, which implies that the spatial 
autocorrelation function can be roughly characterized by C{n) ~ j^-O-'i's 
spatial scales in the range from one to 200 cells. Further numerical analysis 
of the same model setup indicates that if larger spatial scales are considered, 
the power spectrum displays a crossover to white noise at low frequencies. 
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Figure 1: (a): Log-log plot of the power spectrum of a temperature signal 
generated by the term 1/3 (t) of model The parameters are as discussed 
in the text. The power spectrum is consistent with the power-law S{f) ~ 
j-o.25_ xhis indicates a power-law decay of the autocorrelation function 
given by C{t) ~ for time scales between five and 1000 time units, 

(b): An alternative way to estimate the autocorrelation function for the 
same situation as in (a). The standard deviation as a function of time of 
the sum of the temperature signal is plotted in a log-log graph. The graph 
shows a clear difference with respect to the behavior expected from a random 
walk. The standard deviation is consistent with a{t) ~ t~^'^^, which indicates 
C{t) ~ for time scales up to 1000 time units. 
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Figure 2: Log-log plot of the power spectrum of a temperature signal gener- 
ated by the term X3{t) of model (j21). The parameters are the same as in Fig. 
1, except for d2 = 0. A power-law S{f) oc /^°'^ is plotted for comparision. 
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Figure 3: Log-log plot of the power spectrum of the spatial vector of tem- 
perature values vector at a fixed time. The parameters are as discussed in 
the text. For this situation, in which there is no local variability in the pa- 
rameters (besides 7^), the power spectrum indicates a power-law decay of 
the spatial autocorrelation function for scales in the range from one to 200 
cells, characterized by C{n) ~ 

In a more realistic description, the parameters (besides 7„) must have 
some local variability due, for instance, to differences in the Earth's albedo 
and in the solar radiation flux over different regions. As an example, a case 
in which An = A + en and Bn = B + en is discussed in Fig. Hl^a). The e's 
are independent random variables uniformly distributed in the range (—1,1) 
W/m^. The other parameters are chosen as before. The power spectrum 
presented in Fig. Ufa) indicates a crossover between two different scaling 
regimes. For scales from 50 to 1000 cells the correlation function is consis- 
tent with C{n) ~ ?7,~°'^^, while for shorter scales C{n) ~ n~°'^. A situation in 
which the coupling parameters have also local variability is presented in Fig. 
m^b). The values of 6„ are set as 6„ = 6+p„, where p„ is uniformly distributed 
over the interval (—0.15,0.15). The parameters taken as a„ = 0.98 

for all n, which implies a stronger mean coupling. The other parameters are 
the same as in Fig. El^a). In Fig. El^b) the correlation decay is faster for the 
short scales and slower for the large scales in comparision to Fig. Ufa). This 
effect can be interpreted as an increment of large scale coherence as the mean 
coupling grows, while the spatial coherence at short scales decreases due to 
the increment in the local variability. The scaling of the spatial autocorre- 
lation function displayed by the model is reminiscent of recent observations 
on the spatial structure of the climate network, which indicate that nodes in 
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the climate system conform to a network with the small-world property ^H] • 
This property is related to the presence of significant correlations between 
distant nodes. 



(a) 



(b) 





Figure 4: This situation is similar to Fig. except from the fact that some 
parameters besides the 7n's have local variability. In (b) the local variability 
and the mean coupling parameter a have larger values than in (a). 



The numerical findings strongly suggest that model Q shows qualitative 
features that are close to the observations. However, a more precise definition 
of the parameters is needed. For instance, a different choice of the scale 
separation between the 7's and A in general lead to different properties of the 
autocorrelations. Another aspect to be refined is concerned with the already 
mentioned spatial variability of the parameters. This question is closely 
related to the definition of the size associated to cells. In a realistic model 
setup, the parameter values come from spatial averages over the region n. In 
the simple ld+1 situation discussed here, those would be global averages over 
a given latitudinal interval. As already mentioned, these and other relevant 
questions are intended to be investigated by the author in the near future, 
working in close contact with climate experts. 

In order to gain insight into model (j21) it is discussed analytically its 
reduced version given by Eq. 0. Without loss of generality, the constants 
A and B are chosen equal to zero. The term e{t) is a Gaussian white noise, 
defined through the moments {e{t)) = 0, (^e{t)e{t')'j = D6{t — t') and with 
all higher moments equal to zero. The diffusion constant is a parameter 
that measures the strength of the noise. The function 6{t — t) is a Dirac's 
delta. In the absence of coupling x and y simply converge exponentially to 
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the stationary state {y) = x = 0, (y^) = ^ with characteristic times j and 
i. It is assumed that A > 7 so one of the dynamics is fast and the other is 
slow. Strictly speaking, the system has three time scales, the third one being 
associated with the noise. However this time scale has an infinite separation 
with respect to the other two. In the language of control theory, we can view 
y like an output system with uncertainties (noise) that has feedback with an 
input whose response time is different from the time scale of y. Applying 
a Laplace transform over time to Eq. Q and assuming for simplicity the 
initial conditions x(0) = y{0) = 0, one gets 

sy{s) = —\y{s) + \ax{s) + e{s), (5) 
sx{s) = — 7x(s) +'yby{s). 

Solving Eq. © for y{s) and x{s), and by the use of the Faltung theorem, a 
formal solution for y{t) in terms of the noise is found. This solution can be 
used to write explicitly the covariance function, that describes the fluctua- 
tions around the mean value. ^ For large times, the covariance function is 
given by the following expression: 

{y{t)y{t + T)) = q^e'^^^ + q,e'^'^, (6) 

where 

^ /X o,A/^^-7 + 2i^ 'y-X + 2K\ _ 



1^1 = -^^-K, (9) 



/i2 = + K, (10) 



"'^This covariance function can also be derived from a general initial value problem for 
a linear system of stochastic differential equations of arbitrary dimension. 
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K 



\ 



'l + A' 



7A(1 -a6). 



(11) 



In the region of interest of the parameter space, correlations decay monoton- 
ically with a characteristic time 



1 



4^)-\/(4^)"-7A(l-a6) 



7+A 



(12) 



It must be realized that r is greater than any of the two intrinsic times. 
Therefore, because of the feedback there is an emergence of large memory. 
The characteristic time becomes infinitely large when 7 — > with A finite. 
This result is consistent with previous works on reduced models of tempera- 
ture fluctuations 
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Figure 5: Log-log plot of the covariance function Eq. (jHl) with parameters as 
discussed in the text. The covariance cannot be fitted to a single exponential 
for a time interval greater than any of the two intrinsic times. 

The reduced model typically shows a region in which none of the two 
intrinsic time scales is dominant and correlations cannot be adequately fitted 
by a single exponential. For instance, with 7 = 0.2, A = 1, a = 0.96, h = 0.64 
and D = 1(-^)^, the covariance function is not exponential for time scales 
approximately an order of magnitude greater than the intrinsic time of the 
fast variable, as Fig. El shows. Approximation to a power-law or other types 
of slow decay by a sum of exponentials with different characteristic times 
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has been discussed in several fields |19| l2Uj. In particular, this mechanism 
has been already proposed in [T3] , in order to explain the persistence found 
in the atmosphere's temperature record, by fitting the coefficients of a 3d 
AR(1) type process to data. 

3 Conclusions 

The features shown by the spatially extended stochastic model presented 
here motivates the construction of realistic and simple algorithms for the 
prediction of the Earth's temperature distribution and fluctuations. In this 
spatially extended situation the parameters vary locally, so there is a number 
of characteristic times. In order to gain insight on the extended model, 
a reduced version of it has been constructed and the covariance function 
explicitly written. 

One of the future directions of the work is to conduct a more general 
study of the spatially extended model, in close connection to climate research 
to have plausible parameter values. A study of the presented model in the 
context of general systems with several time scales is also intended. The study 
of such systems is important in fields like control theory, inhomogeneous 
media and predator-pray systems among others [TJ |21 El IH El • 
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